To be submitted for publication, 2002 



Transparency of Magnetized Plasma at Cyclotron Frequency 

G. Shvets 

Princeton Plasma Physics Laboratory, Princeton University, Princeton, NJ 08543 

J. S. Wurtele 

University of California, Berkeley, CA 94720 

Electromagnetic radiation is strongly absorbed by the magnetized plasma if its frequency equals 
the cyclotron frequency of plasma electrons. It is demonstrated that absorption can be completely 
canceled in the presence of a second radiation beam, or even a magnetostatic field of an undula- 
tor, resulting in plasma transparency at the cyclotron frequency. This effect is reminiscent of the 
electromagnetically-induced transparency (EIT) of the three-level atomic systems, except that it oc- 
curs in a completely classical plasma. Also, because of the complexity of the classical plasma, index 
of refraction at cyclotron frequency differs from unity. Potential applications of the EIT in plasma 
include selective plasma heating, electromagnetic control of the index of refraction, and electron/ion 
acceleration. 



Electromagnetically induced transparency (EIT) in 
quantum-mechanical atomic systems is a well under- 
stood and thoroughly studied subject. EIT is the 
basis of several very important applications, such as 
slow light ||, information transfer between matter and 
light H,[|, sound wave generation ||, or even testing 
of the black-hole physics Q . Several recent reviews ||] 
illucidated the quantum mechanical mechanism of EIT 
which relics on the destructive interference between sev- 
eral pathways which connect the ground and excited 
states of the atom. The purpose of this Letter is to de- 
scribe EIT in a classical plasma. 

We consider an externally magnetized plasma with 
B = Boe z and density uq. A right-hand polarized elec- 
tromagnetic wave (which we refer to as the probe) at the 
frequency uj\ equal to cyclotron frequency Qo = eBo /mc 
cannot propagate in the plasma because it undergoes 
resonant cyclotron absorption [pi. The cold magne- 
tized plasma dispersion relation ll>i v. s.fci for the right- 
hand polarized probe, plotted in Fig. [H is given by 
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lo 1 = k l c H . Plasma current and the wavenum- 

ber ki become infinite for uj\ — > S7o, and a forbidden 

bandgap develops between Q Q and fi c = fig/4 + w 2 + 

ilo/2, where uj p = (47re 2 no/m) 1 / 2 is the plasma fre- 
quency. This Letter demonstrates that by adding a 
second intense electromagnetic wave (pump) with fre- 
quency luq — CIq — Lo p create a transparency near the 
cyclotron frequency. Moreover, if uj p = Q c , transparency 
can be created by a magnetostatic undulator with arbi- 
trary wavenumber kg. 

The classical mechanism of the electromagnetically in- 
duced transparency is the destructive interference be- 
tween the electric field of the probe Ei± and the side- 
bands of the electric Eq± and magnetic Bq± fields of 
the pump which are produced by the collective electron 
plasma oscillation with frequency lu p along the magnetic 
field. Qualitatively, the total force at the cyclotron fre- 



quency experienced by a plasma electron is given by 
Fiat ~ -e(Ei± + Czd z E ± + C 2 e z X B ±, where ( z is 
the electron displacement in the plasma wave. If the 
pump, probe, and plasma waves are properly phased, 
then Therefore, if the amplitudes and phases of the 
pump and the plasma wave are properly correlated, then 
-Ftot = 0. Consequently, the plasma current at the cy- 
clotron frequency is small (or even vanishing), and the 
probe propagates as if in vacuum. Our numerical simula- 
tion below demonstrates that this correlation is naturally 
achieved in a collisionless plasma. 
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FIG. 1. Dispersion curve for a right-hand polarized wave 
propagating along magnetic field. Forbiden gap exists be- 
tween cyclotron frequency flo = eBo /mc and cutoff frequency 

n c = Q /2 + y/nyi + ui 

We assume two right-hand polarized EM waves prop- 
agating along z— direction, with their electric and mag- 
netic fields given by 2eE ±/mcuJo = a pump e + exp (iSo) + 
c. a, 2eEi±/mcu!i = a pro b e e*+ exp (i0\) + c. c, and 

B 0A± = (cfc ,iM),i) x -^0,1, where e_ + = e x + ie y , 
e_ = e x — ie y , 0q = kgz — ujQt, and 6± = k\z — uj\t. 
Non-relativistic equation of motion of a plasma electron 
in the combined field is given by 
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-j-^ + Q v x e z + ujp( z e z = 



- E 



v x B± r 



m=0,l 



(1) 



where x = {z + ( z )e z + x+ and v — dx/dt are the particle 
position and velocity, and the initial conditions are v = 
and x — z^e z . The third term in the lhs of Eq. (|l|) is the 
restoring force of the ions || . 

Equation (Q) was integrated for two cases: (a) when 
only a probe field is turned on, and (b,c) when both the 
pump and the probe are turned on. The pump and the 
probe amplitudes were increased adiabatically in time, 
up to their respective peak amplitudes of ao and a\, ac- 
cording to 



a 



flprobc = y (1 + tanh [(fi„* - 320)/40]) : 



1 + tanh p t - 160)/40]), 



(2) 



enabling the pump to turn on first, followed by the probe. 

Simulation results for lu p /CIo — 0.3 (u>q = 0.7Qo) are 
shown in Fig. |[ Without the pump electron is resonantly 
driven by the probe as shown in Fig. ||(a). In the plasma, 
this growth manifests itself in a large electron current and 
probe absorption because, time-averaged, E± ■ v± < 0. 
Adding a strong pump with ao = 0.1 and fco ~ 0.83i7o/c 
dramatically changes electron motion, see Fig. ||(b). Af- 
ter the pump is turned on but before the turning on of 
the probe, electron oscillates in the field of the pump 
according to /3 x0 = uj a pump / (wo - Q ) sin (k z - w t). 
Switching on the probe does not significantly alter elec- 
tron motion: (3 X — (3 X $ is shown as a barely visible dashed 
line in Fig. ||(b)). Comparing Figs, ^(a) and (b), observe 
that the pump suppressed electron responce at the cy- 
clotron frequency, making the plasma transparent to the 
probe. 
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FIG. 2. Numerical simulation of the single particle 
motion in the combined field of two EM waves with 
(uji = Qo-M — ioi/c) and (uo = f2o — u P , ko ~ 0.83f2o/c). 
Both pump and probe are slowly turned on according to 
Eq. (|2j). (a) Without the pump electron is resonantly driven 
by probe: (5 X growth indefinitely; (b) With the pump, elec- 
tron motion is almost unaffected by the probe. Solid line 
- total f3 x ; barely visible dashed line - (f3 x — /3 x o), where 
f3 x0 = Ldoa P ump/(uo - Slo)sin(fc 2o — tuot). Since [3 X — f3 x0 . 
(c) Solid line: longitudinal displacement £lo( z /c; dashed line: 
fio«z - Co)/c, where Co = 2a pum p/a probc sin oj p t from Eq. (|) 

Suppression is caused by the excitation of a strong 
plasma [shown in Fig. ^(c)] which produced a sideband 
of the pump at the cyclotron frequency. This sideband 
canceled the electric field of the probe. An approximate 
analytic formula for the steady-state amplitude of the 
plasma oscillation, 



2<2p ro ] :)C . 

Co = t sincjpt, 

^O^pump 



(3) 



is derived below by requiring that the sideband cancels 
the probe, which is in good agreement with the simu- 
lation result. Simulation results demonstrate stability of 
the steady-state values of f3 x and Cz which are naturally 
reached in a collisionlcss plasma. Note that the pump 
has to be switched on prior to the arrival of the probe. 
In atomic physics, this pulse sequence is referred to as 
"counter- intuitive" Q . 

Maintaining high-power pumping waves in the plasma 
may prove challenging in practice. For example, sup- 
porting ao = 0.01 over an area A = (27rc/wo) 2 re- 
quires microwave power of 3 megawatts. Fortunately, 
for Lu p = Qq, a magnetostatic helical undulator can re- 
place a microwave beam. We simulated electron motion 
in the combined field of an undulator with ao = 0.1 and 
ko = 2ilo/c, and a probe which is switched on accord- 
ing to a pro bc = 0.5ai (1 + tanh {(flat — 270)/60]), where 
ai = 0.01. Suppression of the electron responce at the 
cyclotron frequency is apparent from Fig. |§|(a). Electric 
field of the probe is canceled by the (( z /c)e z x B ± force 
which is exerted on a longitudinal plasma wave by the 
helical magnetic field of the undulator. 

Steady-state values of [3 + — (}*_= (3 X — i[3y and Q z 
can be analytically obtained by linearizing Eq. (fil) in the 
weak probe ai ao limit. 

p+ + m p + = 

- {ujoa Q e l§ " + uj iai e lSl - k a C z e iSo - k iai Qe i5 ^ . (4) 

Introducing #o.i = fco.i^o — ^o,it and assuming that 
^o,i Cz < 1j exponentials in Eq. (|4|) are expanded as 
e i§ ,i w e^o.i^ + ik 0A Q), yielding 



-c^aie 101 (l + ihQ - fexG/^i) . (5) 
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Longitudinal equation of motion is given by 

dB\ 



Cz + UpCz ~ \v± X B X + CzV± x . 

1 mc \ oz 



where B±(z,t) was expanded as -Bj_(zo + Cz) ~ B±(zq) + 
Czd Zo B±(zo) to first order in £ z . Inserting the expression 
for B±, obtain 



Cz + U% 



(fc a /3- 



o t0 u 



ik 2 C, z f3_a o e i0 ° - ikl^e 101 ) + c. c. 



(6) 



The last term in the RHS of Eq. (Q) will be later dropped 
because it is proportional to the product of two small 
quantities, Cz an d cl\- Note that, unlike the transverse 
velocity j3+ which is excited directly by each of the two 
lasers according to Eq. (0) , plasma waves are excited only 
in the presense of two lasers via the beatwave mecha- 
nisms. 

The physical reason for EIT in plasma is the strong 
coupling between longitudinal and transverse degrees of 
freedom of the plasma electrons. The steady-state solu- 
tion of Eq. ([j]) Cz = 0.5£expz(Afcz — Au)t) + c. c. , where 
Ao; = u>i — ojo and Ak = k\ — ko, is substituted into the 
transverse equation of motion (|^). Retaining the terms 
with exp —iuot and exp —iuj\t dependence results in 



fa 



iw a _ &i9o 



k>0 — ^0 



u>\ — f2o 



ai + ¥a e*. (7) 



Applying Eq. (Q) to the simulated earlier case of iO\ = fig 
and Auj = uj p yields the steady-state amplitude of the 
plasma wave given by Eq. (||) . 

In the general case of u>\ ^ Slo we insert (3 + and /3_ 
into Eq. (jq) yielding 



-Au 2 )C = ic 2 



(<2i + ik (a /2)- 



kiCLi<jjQ 



• HC^O , , 2 

i- ^-|a | 



(8) 



where Q\ — 9q = (fci — k )z a — Acot. Equation (||) is then 
solved for £ which is substituted into Eq. (^) yielding the 
steady-state value of fa : 



fa 



iuj ao 



jo 



UJidi 



^0 — ^0 ^0 

c 2 kluj \a \ 2 {ki/k - 2) + 



■e idl x 



-O 

2{oj 2 - Aoj 2 )(ujo ~ n ) 



c?^|ao| 2 wi + 2(^2 - Aw 2 )(wi - fi ) 



(9) 



where we have neglected terms proportional to the prod- 
uct of laser detuning <5f2 = lj\ — flo from resonance 
and the pump intensity <2q. Qualitatively, the pump in- 
fluence is strong only close to the cyclotron resonance, 



and is negligible far from uji = £1 - From Eq. (||), 
plasma is resonantly driven when the denominator D = 
2(u 2 - Aw 2 )(wi — fl ) + c 2 k 2 \a \ 2 uj 1 vanishes. Close to 
cyclotron resonance D w Alo p (V1 2 r — <5£! 2 ), where £Ir = 
ckoao{Q.o/ Auip) 1 / 2 is the effective Rabi frequency. Hence, 
the modified plasma resonances are shifted from oj\ = Qq 
to U\ = Qq ± 
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FIG. 3. Same as Fig. 2, except u v — flo, ujo = 0, 
ko ~ 2Q.o /c (static helical undulator is switched on from 
the start), (a) Transverse velocity f3 x and (b) longitudinal 
displacement f2oCz/ c during and after the turning on of the 
probe. 



Fluid velocity component b + w (3^ 



fa 



ikiCzfa 



proportional to exp id\ 



dz(Czfas) 

is given 



(5n + Sn )/(n 2 R -Sfl 2 ) where <5fi (fci) 



by 



(2fl 2 R u; /uj p flo)(k 1 /ko - I). 

Dispersion relation for classical EIT in magnetized 
plasma is derived from the wave equation for the 
probe — (c 2 d 2 — d 2 )E = AndtJ, where the rhs is equal to 
— 4:TricLUienob+ + c. c: 



2 / 2 

z k 



£Q + tffioOi) 



* - w p Wi n 2 R (sny 



(10) 



where it was assumed that the frequency of the pump 
is fixed at lu = Q — lu p . Complete transparency 
(lui = k\c 2 ) is achieved at uii = fl — <5f2o? where 
5D,q rs (2uioH, R /u> p ilo)(rio/koc — 1). Note that this fre- 
quency shift is in general very small in the most interest- 
ing regime of flu <§; uj p : \Sflo\ < AVL 2 r /uj p <C £Ir, and can 
be even smaller near cyclotron resonance when pump and 
probe co-propagate. Equation (^) reduces to the disper- 
sion relation for a single probe in magnetized plasma for 
large detunings (5£l) 2 ^> fl R . The influence of the pump 
is significant only in the vicinity of <50 = 0. 

Note that index of refraction is not identically equal 
to unity at the cyclotron resonance. This is different 
from the quantum-mechanical result for a three-level sys- 
tem pd| , where uo\ — k\c on resonance. It can be demon- 
strated that this difference occurs because multiple Lan- 
dau levels E n = nhflo and corresponding Raman-shifted 
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levels En n — E n + huj p participate in the classical EIT. 
Dispersion relation given by Eq. (10) is plotted in 



Fig. (0) for the same plasma parameters as in Fig. (|l]), 
plus a co-propagating pump with — 0.5lo p . The 
flat band between the Qq ± O resonant frequencies is 
a novel feature which is not present without the pump 
(compare with Fig. (|l|)). The width of this EIT band 
proportional to Qr oc ao can become very narrow for 
low pump amplitude. The corresponding "group veloc- 
ity" (understood in a stricly geometrical sense explained 
below) v g — duoi/dki w 2cil 2 R /ujp can also be made arbi- 
trarily small. Slowly propagating wavepacket of electro- 
magnetic waves is a classical analog of the "slow light" 
in atomic systems ||. 

Qualitatively, the spectacular slowing down of EM 
waves in the EIT plasma can be understood by con- 
sidering the entrance of a probe beam of duration Lq 
into the plasma. In steady state inside the plasma, 
the "slow light" wavepacket of length L f consists of the 
transversely polarized field of the probe \E\\ = \B\\ = 
aimcuji/e and the longitudinal electric field of the plasma 
wave E z = 47reno(2ai/fco a o)- As the pulse enters the 
plasma, it loses photons to the pump at the same rate 
as new plasmons are created (according to the Manley- 
Rowe relation). Classical photon density of a field with 
frequency u> is proportional to the action density oc U/u, 
where U is the energy density. We calculate that the 
ratio of the plasmon to photon density inside the "slow 
light" pulse, 

,2 



UpUs/tOp _ £Iq E z 



20| : 



is 3> 1 if flfj <C ujp. Thus, most photons of the original 
pulse are lost to the pump. Since the index of refrac- 
tion remains close to unity, so is the photon energy den- 
sity. Therefore, the loss of photons is due to the spatial 
shortening of the pulse from Lq to Lf — Lq x (2f2^/o;p). 
Because temporal pulse duration does not change, we re- 
cover the previously calculated v g /c = 2£l 2 R /u>^. It is 
precisely in this geometric sense of v g /c — Lf/Lo that 
the group velocity of the slow light is interpreted. v g is 
not related to the speed of individual photons since their 
number is not conserved during the pulse transition into 
the plasma. 



to, = Q_ * !1„ 



en = si — n„ 
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FIG. 4. EIT dispersion curve, £Io/uj p = 3 and 
Qr/u>p = 1/2. Flat band above flo — £Ir up to Q,o + O/?. 
labeled "EIT Band" corresponds to "slow light" and appears 
only in the presence of a pump. 

One interesting application of EIT in magnetized 
plasma is ion acceleration. While laser-plasma accelera- 
tors of electrons Q| have long been considered as a long- 
term alternative to conventional rf cavity-based linacs, 
the field of plasma-based ion accelerators is still in its 
infancy ||13|j , EIT enables one to conceive a short-pulse 
ion accelerator which consists of a "slow light" pulse in 
plasma with approximately equal group and phase veloc- 
ities. Acceleration is accomplished by the longitudinal 
electric field of the plasma wave. Counter-propagating 
geometry is chosen to match the phase and group veloc- 
ities because v p h = <jj p /|fco + k\ « 0.5cuj p /flo. Matching 
v ph = v g yields ao « w p/^o ^ 1- Other types of ac- 
celerators based on the "slow light" which rely on the 
ponderomotive force also appear attractive because the 
ponderomotive force, which scales as the gradient of the 
energy density E z /Lf oc (u! p /flo)Uo/Vg, increases rapidly 
with decreasing group velocity of the probe. 
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